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A method to derive exphcit formulas for an eUiptic 
generahzation of the Jack polynomials 

Edwin Langmann 



Abstract. We review a method providing explicit formulas for the Jack poly- 
nomials. Our method is based on the relation of the Jack polynomials to the 
eigenfunctions of a well-known exactly solvable quantum many-body system 
of Calogero-Sutherland type. We also sketch a generalization of our method 
allowing to find the exact solution of the elliptic generalization of the Calogero- 
Sutherland model. We present the resulting explicit formulas for certain sym- 
metric functions generalizing the Jack polynomials to the elliptic case. 



1. Introduction 

In this paper we explain a method which yields explicit formulas for the Jack 
polynomials; see |St89LlM79| . We found this method by studying integrable quan- 
tum many-body systems of Calogero-Sutherland type fCTl , Su7 j, and our discus- 
sion will be therefore from the quantum integrable systems' point of view. However, 
we made an effort to make this paper also useful to readers interested in symmetric 
polynomials and not so much in the physics interpretation of our results. To explain 
our notation (which is admittedly less elegant than the one used in |M79l ISt89j 
but closer to the one used in physics) we first give a definition of the Jack polynomi- 
als. (The equivalence of our definition and the one in pt89 follows from Theorem 
3.1 in |St89| .) We then recall the relation of the Jack polynomials to the eigen- 
functions of the so-called (quantum) Sutherland model [SuTl , Su72] and explain 
some physics terminology which we use. We conclude this introduction with the 
definition of the elliptic generalization of the Sutherland model |C75L IOP77| and 
an outline of the rest of this paper. 



A definition of the Jack polynomials. The Jack polynomials Jn(z;l/A) 
are symmetric polynomials of > 1 variables z — {zi, Z2, . . . , zn) G labeled by 
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partitions n (i.e. n e with ni > n2 > ■ ■ ■ > un > 0), which are of the form 

N 

(1.1) J„(Z) = ^ «n,mAf™(z), M™(Z)= ^ 

m<n P6Sjvj=l 

with real coefficients Un,m fixed up to some non-zero normaHzation constant t;n,n 
(which we ignore for simphcity), and which can be defined by the property that 
they are eigenfunctions of the differential operator 

1 ^ ^ 92 a; ,2 Q 



(1-2) D^-yz^ 



2A ^ ^ dz^ ^ Zj — Zk dzj 

J — l J j,k=l -* 

for some parameter A > 0; the ordering of partitions which we use is defined as 
follows, m < n if ni + rt2 + . ■ . + tin = mi + m2 + . . . mN and mi + TO2 + . . . rrij < 
rii + n2 + ■ ■ .Tij for all j, and Sn is the permutation group. ^ [Thus our A, N, and 
Mn(z) correspond to 1/a, n and mx{x) in St89 , and we write vectors with N 
components in bold face.] To see that this defines unique symmetric polynomials 
(up to normalization) one can check that 

(1.3) DM„ = KmM^ 

m<n 

for certain real coefficients 6n,m which one can compute, and this implies that one 
can determine the coefficients Wn,m, m < n, recursively from Vn^n and the condition 
DJn = bn.nJn', SBC |St89j . Theorem 3.1 and the discussion thereafter. In princi- 
ple one can compute straightforwardly the Jack polynomials from this definition 
^St89) . however, the recursion relations one thus obtains for the coefficients Wn.m 
are complicated, and it is therefore difficult to solve them explicitly and obtain 
closed formulas. 

The Sutherland model and its relation to the Jack polynomials. The 

Sutherland model is defined by the iV-body Schrodinger operator 

(1.4) H = -j2-^+^ J2 

i=l i l<j<k<N 

where 

(1.5) V{r) 



4sin^(r/2) 



Xj G [— TT, tt], and 

(1.6) 7 = 2A(A-1), A>0. 

This differential operator has a natural physical interpretation as Hamiltonian defin- 
ing a quantum mechanical model of N identical particles moving on a circle of length 
27r and interacting via the two body potentials V. (To be precise, this Hamiltonian 
is the self-adjoint operator on i^([— tt, tt]^) defined by the Friedrichs extension of 
the differential operator H above, and one is only interested in particular eigen- 
functions of H specified below. We set the length of the circle to 2tt only to ease 



^Usually these monomials Afm are defined with a shghtly different normahzation obtained 
sum summing only over the distinct permutation P ISt89l : the choice of normalization is irrelevant 
in our discussion. 
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notation, and an arbitrary length L > of space as in |Su72| could be easily in- 
troduced by rescaling xj {2n/L)xj, H iJ(27r/L)^, etc.) Following pioneering 
work by Calogero solving^ a similar model |C71j . Sutherland found an algorithm 
to compute the eigenvalues and eigenfunctions of H |Su71l ISu72| , and these exact 
solutions have made these models famous among theoretical physicists (since usu- 
ally such quantum models describing interacting particles can only be studied using 
approximation methods; moreover, there are many interesting physics applications 
and variants of these models; see |OP83j for review). 

Sutherland's solution method is closely related to our definition of the Jack 
polynomials above. To be more specific: The well-known ground state of the Suther- 
land model {— eigenfunction of H with the smallest possible eigenvalue) is 

(1.7) *o(x)= n 

^<j<k<N 

where 

(1.8) 6'(r) = sin(r/2), 

and to obtain the other eigenfunctions of Sutherland made the ansatz 

(1.9) *(x) = $(x)*o(x) 

transforming the eigenvalue equation H'^ — £'i> to H'^ — f $ where H' — 
[l/*o(x)]i?*o(x) equals 

(1.10) H ^'^ ^ gia;^ _ gixfc [qj;, dxk 

j = l 3 j<k ^ ■> 

and £' = £ - Eq with Eq = X^N^N^ - 1)/12 the ground state energy defined 
through -ff^o = -£'0^0- Sutherland's method to compute the eigenfunctions $ of 
H' is equivalent to diagonalizing the operator D in H1.2|l by making the ansatz in 
(|1.1|) etc., as discussed above. To see this, we change variables to Zj = e'^^ and 
obtain 

N Q 

(1.11) F'(x)=2AD(z) + [l-A(Ar-l)]P(z), P^J^z^ — . 

j=i 3 

The operators P and D commute, and it is not difficult to see that the eigenfunc- 
tions $n(x) = Jn(z; 1/A) obtained by Sutherland's method are equal to the Jack 
polynomials. We thus can conclude that the eigenfunctions of the Sutherland model 
are given by 

(1.12) vl/(x) = e'PS"i^^J„(z;l/A)*o(x), = e'^^ 

with arbitrary p G M (wc used the fact that, if ^'(x) is an eigenfunction of H, then 
exp(ip^^ Xj)^'(x) is an eigenfunction of H as well). The corresponding eigenvalues 

are found to have the following remarkably simple formula, £o{n) = "^jLiln-j +p + 
X{N -I- 1 — 2j)/2]^ jSu72| . It is interesting to note that the exponential factor 
on the r.h.s. in H1.12|l describes the center-of-mass motion of the system, which is 
not very interesting and thus often ignored. However, it should be included if one 
is interested in all eigenfunctions of the model. Moreover, since transformations 



■^By solving a quantum many-body model we mean to determine the eigenvalues and eigen- 
functions of the Hamiltonian defining this model. 
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p —I- p + k and nj nj — k for arbitrary integers k < hn leave the eigenfunction 
invariant, one should restrict p to integers and set un = if one wants to avoid over- 
counting. We also note in passing that the ansatz in (|l.l|l has a natural physical 
interpretation: The monomials M„(z) = X^peSw ®-'^P(^Sj ^Pj^j) ^'^s plane waves 
providing a complete set of eigenfunctions of the non-interacting Hamiltonian H 
with 7 = 0, and the Jack polynomials thus correspond to a linear superposition of 
plane waves. 



The elliptic Calogero-Sutherland model. In this paper we explain a 
method to solve the Sutherland model which, different from Sutherland's method, 
gives fully explicit formulas for the eigenfunctions |L01l IL04dj . We will mainly 
concentrate on the Sutherland model, but a main motivation for writing this paper 
is that our method can be generalized to the well-known elliptic generalization of 
the Sutherland model where the interaction potential in (|1.5f) is replaced by 

(1-13) V^(r) = y ^— — , /3>0 

„^^4sin2[(r + im/3)/2]' 

which is essentially the Weierstrass elliptic p-function with periods 27r and ifi.^ 
This so-called elliptic Calogero-Sutherland (eCS) model is known to be integrable 
|C75l [DP77j . which suggests that it should be possible to solve it. For N = 2 the 
eigenvalue equation of the eCS model is essentially equivalent to the Lame equation 
which was studied extensively at the end of the 19th century; see jWW62] for 
an extensive discussion of the classical results. The problem of finding the general 
solution of the eCS model (without any restrictions on parameters) seems to be 
regarded as open, even though there exist various interesting results in this direction 
|m93l lEl?94l IEFKQSL rFV95a[ IFV95bL 15951 ITOOl IKT02L IfNP03. . As we 
will discuss, our method also provides an explicit solution of the eCS model. In 
the elliptic case we will also obtain solutions as in H1.12|l and (|1.7|) above, only the 
function 6{r) is replaced by 

C30 

(1.14) 6i(r) =sin(r/2) ]J(l-2g2™cos(r) + g^"), g e-^/^ 

m— 1 

which is essentially the Jacobi theta function and the Jn(z) are symmetric 
functions which no longer are polynomials but reduce to the Jack polynomials in 
the limit g i 0. We will sketch how to obtain explicit formulas for these elliptic 
generalization of the Jack polynomials by infinite series.^ 



Plan of the rest of this paper. In the next section we explain our explicit 
solution of the Sutherland model and present our explicit formulas for the Jack 
polynomials. Our arguments are such that they generalize with minor changes to 
the elliptic discussed in Section |3| This section also contains a description 

of our results for the elliptic generalizations of the Jack polynomials. We end with 
a few remarks in Section^ 



^To be precise: V(r) = p(r) + co where cq = (1/12) - (1/2) 52„6n l/*^™li^(/^™-/2) IWW62I . 
4To be precise: e(r) = i9i (r/2)/ [2gi/4 pi^^^ (i _ g2m)j |WW62| . 

'^These q-deformed Jack polynomials are different from tfie Macdonald polynomials |M79| . 
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2. Solution of the Sutherland model 



We start by summarizing our solution in a theorem. The proof of this theorem 
has the character of a derivation which not only proves but, as we hope, also clarifies 
and explains this result. 



Theorem 2.1. For m G Z'\ let 

N 

(2.1) u{z) = n 



c, 27ri^i' 



ni<,<fe<Are(g,-/6) 

nl..=ie(z,/e.)^ 



where 
(2.2) 



e(0 = (i-e), 



and the integration contours Cj are nested circles in the complex plane enclosing 
the unit circle, 

(2.3) Cj : = e^^e'» , -tt < yj < tt 
for some £ > 0. Moreover, for partitions n, let 

(2.4) Pn(z)= J2 an(m)/™(z) 



with 

(2.5) 

where 
(2.6) 



a„(m) = (5(n,m) +^7" ^ H ^^'^ 

s=l ji<ki ui is<fcs i^s 

J(m, n + X;^=i '^rEj,.fcJ 



n^=i ^o(n + ELi '^^Ej.fcJ - £0(11) 



fo(m) =^[m, + iA(7V + l-2j) 



/or a/Z m G Z"'^, Ej^ zs the vector in 1^ with the following components, 

(2.7) (Ejfc)f = ijj — 

for j,kJ=l,2,...,N, and 5{n,m.) -.^ Y[jLi^nj,mj- Then 

(2.8) ^'n(x) = P„(z)*o(x), = e'"^ 

with ^'o(x) defined in j j. 7| ) and is an eigenfunction of the Sutherland Hamil- 

tonian H in ^1-4^ and fl.5p . 

(2.9) i?*„(x) = £n*„(x), 
anrf i/ie corresponding eigenvalue is 

(2.10) 5„ = £o{n). 

It is important to note that the sums in H2.4|l and H2.5|l are finite (i.e. all 
but a finite number of terms in these sums are zero), and the /m are symmetric 
polynomials. From this we can conclude: 

Lemma 2.2. The Pn(z) given in Theorem \2. 1\ are symmetric polynomials. 
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From this one can conclude that the Pn(z) are (essentially; see below) equal 
to the Jack polynomials unless there is an degeneracy, i.e., unless there exists a 
m ^ n with ~ ^j) =0 such that £o{m) = £o{n) |L01| . To resolve this un- 

certainty due to possible degeneracies we have recently checked that the P„ (z) also 
are eigenfunctions of the well-known third order differential operator commuting 
with Sutherland Hamiltonian |OP77| .^ and we thus have convinced ourselves that 
always, 

(2.11) (2l22---2Ar)'=Pn_fce(z) =C„,fcJ„(z), G = (1, 1, . . . , 1) 

with some normalization constant Cn,fc. We thus get an infinite number of explicit 
formulas for each Jack polynomial. We checked in the simplest case N = 2 that the 
fact that the l.h.s. of H2.11|l is independent of k is non-trivial. Thus H2.11|l implies 
an infinite number of non-trivial identities which are essentially the contends of 
Theorem 5.1 in |St89| . 

It is interesting to note that the functions /m given above can be non-zero even 
for certain non-partition m G Z^, and in our proof of Theorem 12.11 it does not 
seem essential to restrict to n's which are partitions. We checked for N = 2 that 
the functions Pn all vanish unless n is a partition, and this is true due to highly 
non-trivial cancellations.'' We do not know how to prove this analytically and for all 
N using only our approach, and the same is true for the question of completeness, 
i.e., whether our construction gives all eigenfunctions or not. However, these facts 
can be proven by comparing our solution with Sutherland's |Su71L [Su72| . 

An important aspect of our method is that we expand our eigenfunctions in 
a set of functions /m which are much more complicated than the monomials Afm- 
However, as we will discuss, there exist fully explicit formulas for the /m. Moreover, 
the /m are much closer to the exact eigenfunctions than the monomials My^ in 
the sense that the coefficients Q!n(m) are much simpler than the coefficients fn,m 
discussed in the introduction, and we therefore can compute their explicit series 
representation in (|2.5() . This later formula is our main result in this paper in 
addition to our previous results in Reference |L01) . 



Proof of Theorem 12. 11 We derive the result stated in Theorem l2.1l in three 
steps. 

Step 1: A remarkable identity. The starting point of our solution method 
is a particular functional identity. 

Lemma 2.3. Let 
(2.12) F(x,y) 



Ui<j<k<N ^i^J - XkYOiyk - ViY 



Yllk=l ^i^j - Vk) 

with 9{r) the function defined in M.t^ and x, y G C"*^. Then 
(2.13) i7(x)P(x;y) =iJ(y)F(x;y), 

where H is the differential operator defined in \1.4\j and \1.5]) acting on different 
arguments x and y, as indicated. 



A 



^E.L., unpublished. 

'^These computations were done with the help of MATHEMATICA. 
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This can be proven by a straightforward but somewhat tedious computation 
using the weU-known identity 

(2.14) cot(a;) cot(y) + cot(x) cot(z) + cot(2/) cot(2;) = 1 if z + y + z = 0; 

see Appendix A in LOl . While this provides an elementary proof, it does not 
explain why this identity is true. We thus mention that we found this identity 
by studying a particular quantum field theory model, and this provides a natural 
physics interpretation of this result; see |L04e| for review. 

Remark 2.4. It is interesting to note that the identity in (|2.13(l can also be 
obtained as a corollary of Proposition 2.1 in Reference |St89| . Moreover, it seems 
that our method is closely related to the solution of the Sutherland model by 
separation of variables KMS03 since (|2.13|l seems to give an alternative proof of 
the crucial Theorem 4.1 in (KMS03] . 



Step 2: Fourier- type transformation of the remarkable identity. It is 

useful to note that the identity in (|2.13l) remains true if we replace i^(x;y) by 

(2.15) i^'(x;y) =ce'^^f=i(^^-^'')F(x;y) 

for arbitrary constants P, c (this is not difficult to verify |L01| V Below we will find 
it convenient to choose 

(2.16) P = -XN/2, c = (2/i)^[^(^-i)/2-^^]. 

The idea now is to trade the variables y for suitable quantum numbers n by 
taking the Fourier transform of H2.13|l (with F replaced by F') with respect to the 
variables y, i.e., apply (27r)~''^ J d^y exp(in • y) with suitable Fourier variables n. 
We need to do this with care since, firstly, the function F'(x; y) and the Hamiltonian 
i?(y) have singularities and branch cuts, and secondly, the function F'(x;y) is not 
periodic in the variables yj but changes by phase factors under yj yj + 27r. As 
we will show below in more detail, the second problem can be accounted for by 
choosing the Fourier modes as 

(2.17) rij = rij + ^X{N + 1 - 2j), Uj £ Z, 

while the first problem is solved by shifting the j/j-integrations in the complex 
plan as follows, yj = tpj + \je with real ipj £ [— tt, tt] and some e > (one can 
take the limit e i 0, but this turns out to be unnecessary). Using the fact that 
l/[4sin^(y/2)] = — Yl^=i^^^'^^ ^^r 3(i/) > 0, straightforward computations leads 
to the following result. 

Lemma 2.5. For all n e Z^, the function 



(2.18) F(x;n) = 



F'{x;y), e>0, 



fij in ^2.17}) , and X G [— tt, tt]^, are well-defined and obey the identities 

(2.19) i/(x)F(x; n) = fo(n)F(x; n) - 7 S,F{x; n + t/Ej^) 

j<k i/SZ 

with £Q{n) and^jk defined in Theorem \2.1\ and 

(2.20) 5*0 = 0, S^^iy and S-^ = Vz^ > 0. 
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Moreover, 

(2.21) /'(x;n) = /„(z)vI/o(x) 

with the functions /„ and defined in Theorem \2.1\ 

Note that the integral in H2.18|l is independent of £ > (due to Cauchy's 
theorem). 

It is not difficult to understand how H2.19|l results from (|2.13|) : the l.h.s. is 
obvious since H{x) commutes with the Fourier transformation, while the r.h.s. 
arises from a simple computation using 

(2.22) i/(y)^_^__^^^5,e-fe-2/0 for 3(y, - yfe) > , 

3=1 ^3 j<k,^ei. 

with the first term in the r.h.s. of (|2.19|) coming from the derivative term and partial 
integrations using — Soin), and the second term comes from the interaction. 

To show that F{x; n) is well-defined we use sin(y/2) = iie"'^/^(l — e'^) to write 
(at this point we make a convenient choice for the constants c, P in H2.15|l l 

(2.23) e--^F'(x; y) = ff ni^,<fc<A. 8(0/6)^ 

with 6(0 = (1 - and = e'^^ where the shifts \{N + 1 - 2j)/2 in the pseudo- 
momenta exactly cancel the terms which would make the integrand non-analytic in 
the ^j. This imphes ifT^ . 

Step 3: Ansatz for eigenfunctions and solution of recursion relation. 
Equation (|2.19() suggests the following ansatz for the eigenfunctions of the Suther- 
land Hamiltonian H, 

(2.24) ^'„(x)= ^ a„(m)F(x;m) 

with the normalization condition 

(2.25) a„(m) (5(m,n) + 0(7). 

This, H2.19|l . and the eigenvalue equation (|2.9|l lead to the following relations for 
the coefficients an{m.), 

(2.26) [fo(m) -£„]a„(m) -7^^5,a„(m-i.E,fe) 7(San)(m) 

j<k v& 

(the last equality defines a convenient shorthand notation), i.e., the latter relations 
imply We now observe that ifT^ and yield 

(2.27) a„(m) = (5(m, n) + 7(i?„a„)(m) 
with the linear operator i?n defined as follows, 

(2.28) (i?„a)(m) \ (§c.)(m); 

F0(m) - tnjn 

here and in the following we use the following convenient notation 

(2.29) , I ^ , := [l-^(m,n)]. ^ 



[fo(m)-£:j„ ■ ' ^ ' [fo(m) - f„] ■ 
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Equation (|2.27(l can be easily solved by iteration setting an^m) = S{m, n),* 

oo oc 

a„(m)=^7^(i?>("))(m) = 5(m,n) + ^7^ J] J2 S,, ■ ■ ■ 

s=0 s=l ji<ki ui£Z 

(2-30) X y y 5. ^(m,n + E;,,..E,„,J ^ 

Setting m = n we obtain Q;n(n) = 1 (since (i?nan)(n) — 0), and this and Equa- 
tions (|2.25|) and 12.26|l for m = n imply the following condition determining the 
eigenvalue 

oo 

£„ - fo(n) - -7(San)(n) = - ^ 7^+' E E ^''^ ' ' ' 

^'■''^ ' .^e. n:.r Po(n + ..eUj - „ ' 

We now recall that S*^ = for < 0, and this simplifies the previous formula 
considerably: it implies that all terms in the sum on the r.h.s. vanish (since the 
Kronecker deltas always give zero), and we thus get the result in H2.10|l . Moreover, 
for the same reason we can replace all double brackets |- • - In in H2.30|l by normal 
bracket [• • • ] (since we only have the expressions in (|2.29l) with m 7^ n) . Thus the 
series for Q;n(ni) simplifies to the expression given in H2.5|l . □ 

A simpler argument to derive (|2.10() is as follows: since S^kq = 0, H2.26|l has an 
obvious triangular structure 'IjOI,, and inserting m = n thus immediately implies 
(|2.1U|I . However, this argument does not generalize to the elliptic case, whereas 
ours does. 

Proof of Lemma 12.21 The functions /„ in H2.1|l and (|2.2(l can be computed 
by Taylor expanding the integrand using the binomial series and computing the 
^j-integrals which project out the -independent terms of the integrand. One thus 
obtains (for details see |L01| . Appendix B.3), 

(2.32) /„(z) -ypn,„,M™(z) 

m 

with the Mm(z) in and the coefficients are 

(2.3.) n CJn 

where the sum ^ here is over all non-negative integers fJ.j'k',t^jk restricted by the 
following 2A'' equations, 

N j-1 N JV 

(2.34) uj = y i^ij + y ^lij - y = y 

1=1 1=1 e=j+i e=i 
and TOi > TO2 > . . . > TTiAT > 0. This shows that the /„ are symmetric polynomials 
which are non-zero only if 

(2.35) nj+nj+i + ... + nN>0 Vj = 1, 2, . . . TV, 

■^We realized this only when rereading our paper .LOl. during this spring. 
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and the sum in H2.34|l only contains terms with 

N N 

(2.36) ^mj=^nj. 

i=i j=i 

The latter implies that all series in H2.32|l and H2.5|l are finite (i.e., they truncate 
after a finite number of terms). We finally need to show that the denominators in 
(|2.5|l are always non-zero. For that we write 

(2.37) m = n + fi with fi = fijkEjk, 

and with that and H2.6|) we compute 
(2.38) 

fo(m) - £:o(n) = XI M H ^J^t^i - nt + (k - j)X] + ^ikJ - fJ-jk 

j=i y k=j+i k<j k>j 

which is manifestly positive if n is a partition and all fXjk > with fi ^ Q. Since 
obviously only such terms appear in the denominators in (|2.30() we conclude that 
all Q;n(m) in H2.5|l are well-defined finite series. □ 

As already mentioned, in our proof it does not seem essential to restrict to 
n's which are partitions. In fact, we only used this property to have all energy 
differences £o{ra) — i?o(n) in H2.38|l manifestly positive, but we only would have 
needed that these energy differences are always non-zero which is true also for 
certain non-partition n's, e.g., for N = 2 and A non-integer. We feel that this point 
would deserve a better understanding. 

3. Solution of the elliptic Calogero-Sutherland model: Outline 

We formulated Theorem 12.11 and its proof so that it straightforwardly extends 
to the elliptic case: one only needs to replace V, 9, 8 and 5^ by their elliptic 
generalizations. To be more specific we describe the generalizations of steps 1-3 in 
our proof of Theorem 12. II in more detail. 

Step 1: Lemma \2. yi holds true as it stands with H the eCS Hamiltonian defined 
in ^1-4^ and 1^1.1!^) and the function 9 defined in j j. 

This can be proved by a brute-force computation using the following well-known 
identify for the Weierstrass elliptic functions ^ and p [WW62 , 

(3.1) Ux) + ^{y) + ^{z)f = p{x) + p{y) + p{z) ff x + y + z = 0; 

see Appendix A in |L04b| . A quantum field theory proof of this results can be 
found in |L04a| . 



Step 2: Lemma \2.5\ holds true as it stands with Si, in \2.2(]\) replaced by 

1 

1 — (7^'^ 1 — 

and the functions /„ and ^Pq in 1^2.1]) and jJ.7| ) with 



(3.2) 5*0 = 0, Sy = V- — and S-i, — v- ^ Vi^ > 



(3.3) 9(0 = (1 - n (1 - - 9""/^) 
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and in jj.l^l ). 

The proof is essentially unchanged, only now the functions 8 and the coeffi- 
cients 5*1, in the identities 

(3.4) 9{y) ^ iie-'^/2Q(giy) 
and 

oo 

(3.5) V{y) S.^'^y for 3(zj) > 

need to be changed to what is given above, and we also need to assume e < (3/N; 
see Appendix B in |L01| for a detailed proof. 

Step 3: The argument given remains unchanged until and including l|2.31|l . 
The crucial difference now is that Si, no longer vanishes for negative values of v, 
and thus the above-mentioned triangular structure is lost. Due to this (|2.31(l does 
not simplify but remains as an implicit equation determining the eigenvalue £n, 
and thus the eigenvalues for q > are much more complicated. 

We can summarize the solution of the eCS model thus obtained as follows. 

Theorem 3.1. The eigenvalues of the eCS model are determined by 
with S,y in and £"0(111) in 1^2. The corresponding eigenf unctions are given 

m with F (x; n) in J^tjDjl . {T^j, and (TT^ , and the coefficients 

an(m) m WHIj) . 

It is interesting to note that H2.31(l can be turned into a fully explicit formula 
for the eigenvalues as follows: Defining the function 



(3.6) X 

rir^i [[^o(n + ELi vi^i^k,) - C 

of one complex variable ^, we can write H2.31|l as 

(3.7) f„ = £o(n) + $n(f„). 

Using Lagrange's theorem as stated in |WW62] . Paragraph 7.32, the latter equa- 
tion can be solved by the following infinite series, 

(3.8) £:„ = f„(„)+^_ <i>„(^)^ 

"=1 C=£o(n)+a 

where a is a real parameter which formally can be set to any value such that <&n(C) 
is non-singular in ^ = £o{n) + a. From the explicit formula for <i>n(C) given above 
it is straightforward to compute the coefficients in this series explicitly |L04d| . In 
a similar manner one can find fully explicit formulas for all coefficients an(in). 

The arguments sketched above are enough to obtain an explicit solution in the 
sense of formal power series in q. However, it is obviously important to also check 
if the analyticity properties of the function $n(0 such that Lagranges' theorem 
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gives convergent series. For that the parameter a in H3.8|l is important: One can 
prove that, if a is such that 

N 

(3.9) 3A > : Vm 7^ n with ^(toj - n^) = : \£o{m) - £o{n) - a\ > A, 

then the function <i>n(0 analytic in a disc around ^ = ^o(n) + a which, in a finite 
g-interval, is large enough for the series in (|3.8|) to converge. A similar results holds 
true for all the coefficients an(m). We are confident now to be able to prove that 
the eigenfunctions are square integrable in a finite g-interval as well. 

The parameter A is obviously important since it restricts the range of conver- 
gence of our series solution. It is not easy to give general lower bounds for A, and 
we only mention one simple case: if A is integer then A > 1/2 for a = 1/2. It is 
possible to increase A by finding 'better' values of a, and we believe that this is 
useful for computing numerically using our results. 

We plan to give a more detailed discussion of this solution of the eCS model 
together with detailed proofs in a future revision of Reference jL04b| . 



4. Final remarks 

1. It would be desirable to compare our perturbative solution of the eCS model 
with brute-force numeric solutions and thus check if our results arc also numerically 
useful.^ 



2. As discussed, the eigenvalues of the eCS model are determined by the 
implicit equation in H3.7|l . While the explicit series solution in (|3.8|) can be regarded 
as perturbative solutions which continuously deforms the solution for g = 0, we 
do not see any reason to rule out the possibility that there are additional non- 
perturbative solutions for larger values of q. The spectrum of the eCS Hamiltonian 
as a function of q might therefore be quite complicated with qualitative changes at 
certain critical values of q (this should be closely related to the intriguing analytical 
structure of the functions $n('^) defined above). If so the eCS model would challenge 
notions of quantum integrability based on the simplicity of the spectrum. 

3. Lemma 12.21 has an interesting generalization. 
Lemma 4.1. For non-negative integers N,M, let 

ni<j<fe<Ar ^{Xj - Xk)^ Ul<j<k<M HVk - Vj)^ 



(4.1) FMM{^,y) = 



en 



with 9{r) %n and x e C^, y £ C*^ Th. 

(4.2) [Hx,Ni^) - i^A.M(y) - cn,m] FN,Mi^-. y) = 

where H = Hx,n{^) is the differential operator defined in \1.4\j - J^l.b]) and Cpf^M = 
X^{N - M)[{N - Mf - 1]/12. Similarly, the function 

N M 

(4.3) fw,m(x, y) = n ^(^^ - ^fc)^ n ^^y^ - w)'^' n n ^(^^ - 

l<j<k<N l<j<k<M i = lk=l 



thank David Gomez-UUate for stressing this point. 
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obeys the identity 

(4.4) [i?A,7v(x) + \Hi,x.M{y) - cn.m] Fn.m i^;y) = 

with CNM = [X^NiN"^ - 1) + M(M2 - 1)/A + 3MN{XN + M)]/12. 

We thus can have different particle numbers for the two sets of variables x and 
y, and in addition to the case where the coupling A is the same there is also a 
dual relation with reciprocal couplings for the x and y variables. Note that, for 
AI = 0, H4.2|l and (|4.4|l both reduce to the eigenvalue equation H'^o = Eq'^o for 
the groundstate in H1.7|l . 

We found these identities using a quantum field theory construction, but once 
known they can also be proven by brute- force computations ,LQ4e. Taking the 
Fourier transform of 14. 2|) and (|4.4(l etc. (as in Section|21), one can obtain additional 
explicit formulas, and thus identities, for the Jack polynomials. We suspect that 
(|4.4|) would lead to an alternative proof of Theorem 3.3 (duality relation of the Jack 
polynomials) in |St89| . 

It is important to note that the generalizations of (|4.2|) and (|4.4|) to the elliptic 
case exist but involve, in general, a term with a /3-derivative L04c , and this term 
is only absent in the special case stated in Lemma 12.31 We only know how to 
construct eigenfunctions when this /3-derivative term is absent, which is the reason 
why we only discussed this case in detail. 

4. The method explained in this paper can be adapted to give explicit solu- 
tions of the original Calogero model |C71j and the BCAr-variants of the Calogero- 
Suthcrland models OP83 It will be interesting to see if this method can also 
provide an alternative solution of the Ruijsenaars models iR87j and thus can give 
interesting explicit formulas for Macdonald's generalization of the Jack polynomials 
|M79j : see e.g. |AOS96| for interesting results in this direction. 
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